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CALCULATION OF CURRENT COLLECTED IN A DILUTE PLASMA THROUGH 
A PINHOLE IN THE INSULATION COVERING A HIGH-VOLTAGE SURFACE 

by Gustave C. Fra lick 
Lewis Research Center 

SUMMARY 

A procedure is described for calculating the current collected through a pinhole 
defect in the insulation covering a high-voltage surface. The results apply to a satellite 
at geosynchronous altitude where collective plasma effects and the motion of the satel- 
lite can be neglected in computing the currents. 

The results indicate that the presence of the insulation around a pinhole reduces the 
current and that the experimentally observed current enhancement is probably due to 
surface leakage in the insulation. 


INTRODUCTION 

The use of high-voltage solar arrays on spacecraft to power devices such as radio- 
frequency amplifiers and ion thrusters promises more efficient spacecraft and less 
complicated electronics. The feasibility of using such arrays has be'en investigated 
(refs. 1 and 2). It was found that the possibility of large power losses due to current 
leakage to the space plasma, especially at lower altitudes, can be minimized by insulat- 
ing the high-voltage surfaces. However, voids in the insulation and puncture by 
micrometeorites will result in high-voltage pinholes being exposed to the space plasma. 
There is evidence that these pinholes can collect unexpectedly large currents. Cole, 
Ogawa, and Sellan (ref. 3) found that the current collected is larger than can be ac- 
counted for by simple analytical models, as have Grier and McKinzie (ref. 4). These 
currents can result in a loss of power and in the formation of a high-temperature region 
at the site of the pinhole, which will lead to further deterioration of the insulation. 
Another effect, not fully understood, is the area effect reported by various observers 
(refs. 5 and 6). The effect of the insulation in these experiments is to actually enhance 
the current collection. The collected current appears to be more a function of the area 



of the insulating surface surrounding the pinhole rather than a function of the pinhole 
area itself. This phenomena has been observed for various types of insulation. 

Most of the analytical work in this area has been done by Parker and Whipple in 
their study of electrostatic probes on satellites (refs. 7 to 9). In their method currents 
are calculated by determining the trajectories of collected particles. They used two 
models for their calculations of probe currents. Their first model was a small circular 
equipotential surface, which represented the probe, surrounded by an infinite flat sur- 
face at zero potential, which represented the remainder of the satellite. In their second 
model the satellite was represented by a right circular cylinder with the probe being a 
small circular area in the center of one end of the cylinder. Published results using 
this model consider only the case where each surface is an equipotential. 

In this report the effects of the insulation are incorporated by modeling the insulator 
with a nonequipotential surface surrounding the pinhole, which should be more realistic 
than using an equipotential surface. In addition, higher pinhole voltages are used, up to 
10 kilovolts, rather than the lower voltages previously considered. These voltages are 
more representative of actual values. The calculations of the currents were performed 
using the trajectory tracing technique of Parker and Whipple. The goal of this study is 
to determine the effects of a nonequipotential surface surrounding the pinhole in the hope 
of providing an explanation of the area effect based solely on electrostatic effects. 


MODEL 

Geometrical Model 

The geometrical model used to simulate a pinhole defect in the insulation covering 
a high voltage surface was simply two concentric disks (fig. 1). The inner disk repre- 
sents the pinhole and is at a uniform potential. The annular region between the radii of 
the inner and outer disk represents the insulation. A typical voltage profile is shown in 
figure 2, which is based on measurements by Boeing (ref. 6). In all the calculations 
the diameter of the outer disk was 10.2 centimeters (4 in. ). Two diameters were used 
for the pinhole, 0. 38 centimeter (0. 15 in. ) and 0. 57 centimeter (0. 225 in. ). This model 
is similar to the pinhole- insulation configurations used in testing in plasma tanks and to 
the test surfaces that were to have been used on the SPHINX (Space Plasma High voltage 
INteraction experiment) satellite (ref. 10). 


Analytical Model 

To calculate the pinhole current, it is necessary to know the voltage on the insula- 
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tion, as this is part of the boundary condition of the problem. But this voltage must be 
calculated. If an insulator is inserted into a plasma, the voltage at each point will 
change until an equal number of positive and negative charges per unit time are attracted 
from the plasma at each and every point on the insulator. If the voltage at a point is 
too positive, an excess number of electrons will be attracted, and the voltage will de- 
crease until the electron and ion current densities are equal. In the same way, if the 
voltage is too negative, more ions will be attracted and the voltage will rise. 

Ideally, then, the current densities should be calculated at each point on the insula- 
tor, and the voltage distribution adjusted until the current densities are equal. The 
voltage on the pinhole would remain constant at its specified value. However, a voltage 
change at one point on the disk would change the voltage distribution everywhere around 
the disk and the current everywhere else, not just at that point. This would lead to 
some sort of iterative technique for computing the insulator potential distribution, which 
may converge only slowly and could be very costly in terms of computer time. 

Instead, it was felt that the magnitude of the effect of insulator voltage could be 
studied without having the "correct” voltage at each point on the insulator. The 
approach was to work with the plots of the equipotential lines surrounding the disk for 
various voltage distributions on the insulator and then to calculate the currents for the 
cases of interest. The plots could be generated relatively quickly. 

The electric field was calculated assuming an infinite Debye length (Laplace field), 
and, although a magnetic field is included in the equations herein, all the currents were 
computed with B = 0. In calculating particle trajectories, it was assumed that colli- 
sions with other particles could be ignored and that the plasma was nonstreaming. Thus, 
the results should be applicable to geosynchronous altitude, where high-voltage solar 
arrays could be used to power communications satellites. At this altitude, the Debye 
length is about 200 centimeters, the mean free path is of the order of several kilometers, 
and the thermal velocities of the ions and electrons are at least an order of magnitude 
greater than the speed of the satellite. 


BRIEF DESCRIPTION OF SOLUTION 

This section contains a short description of how the currents were calculated by 
tracing out particle trajectories and how the electrostatic field was calculated. A full 
description of the method of solution is given in appendix A. (Symbols are defined in 
appendix B. ) 

The normal component of the current density for a given species of particle at a 
point P on the disk is given by 
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j(r p ) = q fff V • nf(r p) V)V 2 dV sin 6 y d 6 y d <p y 
allowed 

where q is the charge on the particle, V is the particle velocity, n is the unit normal 
at the current collection point, f is the distribution function in the six- dimensional 
phase space formed by the three particle coordinates and three components of velocity, 
r p is the position vector of the current collection point, and 6 y and y y give the 
direction of V at P (see appendix A). 

The method of solution makes use of the fact that, although the functional form of f 
is unknown at P, as long as collisions with other particles can be ignored and energy 
is conserved along particle trajectories, the numerical value of f is constant on a 
trajectory (refs. 6 to 8). At infinity, where the presence of the disk has essentially no 
effect on the space plasma, the distribution function is assumed to be a known function 
of the coordinates and velocity. Here, the distribution function is taken to be Maxwell- 
Boltzmann, which, for the case of a nonstreaming plasma considered in this analysis, 
has the form 


*(?«.V eo )=n f 


m 


2nkT 


3/2 -mV“/2kT 
1 e °° =f(r p ,V) 


where is the particle velocity at infinity, n Q is the particle density, m is the 
particle mass, k is Boltzmann's constant, and T is the temperature. The current 
density at P is the result of integrating over only those trajectories that connect with 
the space plasma. To find whether a trajectory is acceptable, a particle is assumed 
to arrive at P with a given kinetic energy and set of impact angles 6 y and <p y . The 
direction of the particle is reversed, and it is sent back away from P, either to end up 
in the space plasma or land elsewhere on the disk. If the particle lands elsewhere on 
the disk, the trajectory is not an allowed trajectory, since no particles are emitted by 
the disk (photoelectric effect and secondary emissions have been neglected). If the 
trajectory connects with the space plasma, it is allowed and contributes to the collected 
current density at P. 

The particle moves along its path under the influence of the electrostatic field and 
of any magnetic field which may be present, according to the equation of motion 

m — = q (E + V x B) 
dt 


The trajectory is found by integrating this equation numerically, using the requirement 


4 



of energy conservation to adjust step size . 

The electrostatic field was found by using the principle of superposition. Closed- 
form solutions exist for the voltage and the components of the electric field around a 
circular disk of uniform voltage. The solutions from several of these uniform voltage 
disks of different sizes are simply added together to get the solution for this problem 
where the voltage on the disk is not uniform (fig. 3). The voltage on each of these disks 
is chosen so that the total voltage from all disks gives the desired voltage distribution 
from the center of the pinhole to the outer edge of the insulation, that is, so that the 
boundary conditions are satisfied. 

In practice, the voltage on the insulator is specified at the inner radius where the 
pinhole and insulation meet, at the outer radius, and at various points in between. A 
typical voltage profile is shown in figure 2, which is based on measurements bv Boeing 
(ref. 6). 

This superposition approach, together with the fact that long-range collective 
effects are not included, avoids the necessity for the finite difference approach used 
by Parker and Whipple (refs. 7 to 9). 

The starting point for calculating the pinhole currents was the voltage profile shown 
in figure 2. This will be referred to as the "very positive" case because of the very 
high positive voltages over its surface. The pinhole had a radius of 0. 19 centimeter 
and was at 10 000 volts. The radius of the dielectric was 5. 1 centimeters. The com- 
puted voltage distribution in the space around the disk is shown in figure 4. 

The electron and ion current densities j g and j. were computed at several points 
on the disk. The results are given in table I for temperatures T = T g = T.. (For high 
altitudes, the ion and electron temperatures are roughly equal. ) The particle temper- 
atures are expressed in terms of E th , the random thermal kinetic energy of the parti- 
cle. If the voltage profile shown in figure 2 had been the correct one for a disk with a 
perfect (no leakage) insulator, then the computed currents j g and j. would have been 
equal at each point on the insulator. 

Instead, the electron current was everywhere too high relative to the ion current, 
which means that the voltage on the insulator was everywhere too positive. Figure 4 
shows that the voltage is quite positive even far from the disk; this case, the "very 
positive case", represents an upper bound on the insulator voltage. 

A discrepancy exists here because although this voltage profile is based on values 
obtained experimentally the calculations indicate that the insulation was too positive. A 
possible explanation lies in the fact that for the purpose of these calculations the insu- 
lation was assumed to have no leakage of current along the surface into the pinhole. As 
opposed to this ideal insulator, however, a real insulator could allow electrons to be 
drawn toward the very positive pinhole, preventing the buildup of negative charge and 
allowing the insulation to remain more positive than it would with no leakage. 
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This leakage could also explain the area effect, since now the total pinhole current 
would consist not only of current collected at the surface of the pinhole, but also the 
current that leaked in along the insulator. The additional current should increase 
roughly with the area of the insulation, which is in fact the case. Further evidence for 
this hypothesis is provided by the fact that an ideal insulator tends to mask the presence 
of the pinhole and reduce the collected current from its value with no surrounding in- 
sulation. 

To test the effect on pinhole current of the voltage on the surrounding insulation, 
three other cases were considered - the slightly positive case, the slightly negative case, 
and the very negative case. These case names are descriptive of the voltages surround- 
ing the disk. The voltage profiles are shown in figure 5, and the voltage contours in 
figures 6 to 8. 

The electron current densities at the center of the pinhole are shown in table n(a) 
for each as (i e /j eo )(E th /> p ), where i// p is the pinhole voltage (10 000 V in all cases). 
The ion current was negligible. 

It is apparent from table n(a) that the current collected by a pinhole is strongly in- 
fluenced by the voltage on the surrounding insulation. This is especially evident in the 
results of the slightly positive and slightly negative cases at low electron temperatures. 

If the presence of the pinhole can be felt even slightly by the electrons, considerable 
current will be collected. But, if the pinhole is masked by the presence of the negative 
insulation, the collected current is greatly reduced. 

The importance of the masking effect of an ideal insulator on collected current was 
also shown in two other ways. 

First, the effect of pinhole size was considered. The pinhole currents were calcu- 
lated for the four cases using a 50 percent larger pinhole. These currents are shown in 
table n(b), and the voltage contours in figures 9 to 12. Comparison of the two positive 
cases for the two pinhole sizes shows that once the pinhole can be "seen”, pinhole size 
does not have much effect on current density. Since the current density was fairly uni- 
form across the pinhole surface, the total pinhole current in this case would be roughly 
proportional to pinhole area. Increasing the pinhole radius in the "slightly negative" 
case, however, resulted in a large increase in the current density. Comparison of th6 
voltage contours in figures 7 and 11 shows that this can be explained by the fact that the 
pinhole voltage "punched through" the masking effect of the surrounding insulation. 

In the second test the current was calculated for a 0. 19- centimeter-radius 10- 
kilovolt circular electrode with no insulation surrounding it. The currents are shown in 
table m, along with the currents for the very positive case. The voltage contours are 
shown in figure 13. Comparing figures 4 and 13 shows that the presence of the insulator 
tends to reduce the voltage in the space surrounding the pinhole, and this reduction in 
voltage is reflected in a reduction in current collected by the pinhole (table III). This is 
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true even in the comparison made here, where the voltage profile used on the disk was 
too positive, as has been shown. Therefore, any area effect must be due to some mech- 
anism not included in these calculations. As discussed artier, the area effect is prob- 
ably due to surface leakage on the surrounding insulation. 

Even if some way is found to reduce this surface leakage, these results indicate 
that pinholes could still be a problem on satellites, where high-voltage surfaces would 
be covered by insulation. A small, isolated pinhole may not have significant leakage to 
the space plasma, but as micrometeorites cause other nearby pinholes, the masking 
effect of the insulation will be reduced, and the leakage current could increase rapidly 
because of the combined effects of increased current density and increased pinhole area. 

For completeness, the effect of different pinhole voltages on the collected current 
was also considered. The results are shown in figure 14. To get the variation of cur- 
rent with voltage, all the voltages on the very positive case were divided by 10 (p -p = 

1 kV) and then by 10 again (ip p = 100 V). The shape of the voltage contours is the same 
as that shown in figure 4, divided by 10 and 100. For j// p >> E th 


j e E th 
j eo 


constant 


Since 


then 


■*eo n e% 



= n. 



1/2 


j e “ (constant ) 1 // p n e 


l3 ” m e E tW 


Hence, if the plasma temperature is much lower than the pinhole voltage, j is directly 
proportional to electron density and pinhole voltage and is inversely proportional to the 
square root of the plasma temperature. The inverse temperature dependence in this 
case is reasonable, as the current collecting surface will have a greater effect on elec- 
trons with less thermal energy, so that the current will increase as the temperature 
decreases. 
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On the other hand, if E th >> i// p , (j e /j e0 )(E th /i// p ) approaches the limit (E th /«// p ), or 
j « j eQ . This is to be expected, as the low voltage will have little effect on the ener- 
getic high- temperature electrons, and the current density will be just the random cur- 
rent density j . Thus, the current calculation techniques used here give reasonable 
results for the variation of current with voltage. 


CONCLUSIONS 

A procedure using the trajectory tracing method of Parker and Whipple has been 
described for calculating the current collected from a dilute plasma through a pinhole 
defect in the insulation covering a high voltage surface. The geometrical model used is 
two concentric circular disks. The inner disk is at a constant potential and represents 
the pinhole. The region between the inner disk and the edge of the outer disk represents 
the surrounding insulation and is not at a uniform voltage. 

Currents are calculated for cases where the voltage profile on the insulation is 
adjusted so the potential field surrounding the pinhole becomes either very positive, 
slightly positive, slightly negative, or very negative. The effect of pinhole size and 
pinhole voltage is also considered. 

The results indicate that the large area of relatively negative insulation tends to 
mask the presence of the pinhole and that both the voltage in the space around the pinhole 
and the current collected by the pinhole are reduced by the presence of the insulation. 
This masking effect means that the ’’area” effect reported by some observers, in which 
the surrounding insulation increases pinhole current, cannot be explained solely by 
electrostatic effects. Instead, the results of this report indicate that the area effect is 
probably due to leakage of current along the insulation surface into the pinhole. 

The masking effect is reduced as pinhole area is increased relative to the insulation 
area. This is likely to occur on satellites as micrometeorites puncture the insulation 
covering high-voltage surfaces. When this happens, all the pinholes could start to col- 
lect large currents. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 5, 1975, 

505-04. 
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APPENDIX A 


METHOD OF SOLUTION 

This appendix includes a detailed description of how the currents are calculated by 
the trajectory tracing technique and how the resulting integrals are evaluated. Also in- 
cluded are descriptions of how the particle trajectories are calculated and how the elec- 
trostatic field around the disk was calculated. 

The normal current density at a point P on the disk whose position vector is F p 
is 


j(r p ) = fff qf(r p ,V)V • nV 2 dV sin 0 y d0 y d<p y (Al) 

allowed 

where 0 y and cpy give the direction of V at P (fig. 15). In equation (Al) ’’allowed" 
refers to the fact that only those particle trajectories that connect with the space plasma 
contribute to the current collected at P. If the current is being collected on the top 
side of the disk, the possible values for 0 y and <Py are 0 ^ (py ^ 2n and 
7t/ 2 < 0 y < it. 

The method for calculating the current makes use of the fact that f is a solution of 
Boltzmann's equation so that it has the same numerical value everywhere along a tra- 
jectory (refs. 7 to 9). Hence, it has the same value at P, at one end of the trajectory, 
as it did at its starting point in the space plasma, at the other end of the trajectory. 

This is true even in a magnetic field, as long as collisions can be ignored so that trajec- 
tories are curves of constant energy. 

The distribution at "infinity" is taken to be Maxwell- Boltzmann and nonstreaming: 


«yv> 


= "o fej 


3/2 -m 

e 


vl/2kT 


(A2) 


It is more convenient to work in terms of energies. Equation (A2) becomes 


\ ,/-» p, , . 3m \ 3/2 _- 3E ke~/ 2E th 

f(r p’ E ke ) ke^ o( 4jr | q | E 


(A3) 
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where 


— mV 2 


rms 


= - kT = | q| E 


th 


(A4) 


~ mV 2 =|,| Ej te 


The random current density in space j Q is the charge crossing a surface per unit 
area per unit time. For a Maxwell- Boltzmann distribution 



(A5) 


Since the collected currents will be written in terms of j Q , equation (A3) is rewrit- 
ten as 


f ( r co > Ev.oo ) - “ ~ 


m 


0O, keoo/ ~ 8 tt q 1 1 q| E 


-3E, /2E., 

keoo' th 


(A6) 


thy 


Making use of the fact that energy is conserved along a trajectory gives 


«r E > 9 J o/m \V 3/2 {K +(q/kl) ' , ' (F P > ] /E th} 
P ’ ke 8, q (|^ e 


(A7) 


where i//(r p ) is the potential (in volts) at P and (q/|q| )^(r ) is the potential energy 
(in eV) a + P. 

It is now possible to use equations (Al) and (A7) to calculate the current collected 
at P. First, however, note that 


and that 


V • n = V • -z = -V cos 0, 


V 3 dV = 


= 2 /IsA 

\m/ 


E ke dE ke 
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It is also necessary to consider the allowed range of integration. For a given value of 
9 y and (p y, there will be a minimum kinetic energy below which no particles will be 
collected. From energy considerations alone, the minimum energy (in eV) of a particle 
collected by an attractive potential must be at least the voltage at the collection point. 

It could be zero for a repulsing potential. The minimum kinetic energy may still be 
higher than that allowed just on the basis of the voltage of the collection point. This 
may be seen by considering a particle ejected from the current collection point. If the 
particle is ejected with too low a kinetic energy, it will fall back somewhere else on the 
disk and not connect with the space plasma (fig. 16). 

At a given point P on the disk, the minimum kinetic energy will depend on fly 
and cpy. For instance, suppose the electron current collected at the center of the pin- 
hole is being calculated and that the potential on the disk is everywhere positive (as in 
fig. 2). The minimum kinetic energy will be greater if 0y is near u/2 than if 0y is 
near ir, because, to avoid being captured at some other point on the disk, an electron 
arriving on a trajectory that just grazes the surface will have to be going faster than one 
coming from straight above. Equation (Al) may now be written as 


i< F p ) = ' J o — T 

4 > E th 


•'<p v = 0 ' 2 ^ 


-3/2 {[ Efce * (q/| q |)«r p )] /feyj 


J ke~ E min( 0 V’ ^ 


x E ke dE ke cos 0 y sin 0 y d0 y dcp y (A8) 
The integral over kinetic energy can be performed analytically by using the result 


f 


-ab 

xe -ax dx = (ab + 1) 


Then equation (A8) becomes 


r n 

io " Jo Jir/2 L 2 


E min( e V’ ‘Py) 


E 


+ 1 


th 


- 3 / 2 {[ E min(V ^yMq/l q| )^( F p)]/ E th} 


X Sin 0y COS 0y d0y d^y (A9) 

The integrations in equation (A9) must be performed numerically. No matter what 
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numerical integration scheme is chosen, it is necessary to know E . for various 

min 

values of 9 y and <Py. The technique used here was weighted Gaussian quadratures, 

r b 

which works in the following way: The integral I f(x) dx is approximated by the sum 

•'a. 

N 

2 WjffXj) 

i=l 

where the Wj are weights assigned to these points and the x i are various values of x 
in the range a < x < b. The values of w t and are found in handbooks (ref. 11) 
and depend on the length of the interval and on N, the number of points used in the ap- 
proximation. Generally, Gaussian quadrature gives the best accuracy with the fewest 
values of the integrand. In this case it meant having to evaluate the fewest values of 
E m in* Each value of E m i n required considerable computer time to evaluate. 

To find E m j n for a given value of 0 y and <Py, consider a particle arriving at 
the current collection point P and having these impact angles. It will have some 
kinetic energy E kg at impact. Since collisions are being ignored, the trajectory is 
dynamically reversible. The particle is sent back along the same path, starting at P 
with E^ e . The path will terminate either in the space plasma or somewhere on the disk. 
If the particle came from the plasma, E kg > E min , so the process is repeated with the 
same 6> y and <p y and a smaller value of E ke until the path terminates on the disk. 
Then E ke s E mkl (see fig. 16). The range between the upper and lower bounds on 
E min is reduced un til a sufficiently accurate value of E m - n is found for this value of 
6y and <p y . The criterion used here was that the difference between the upper and 
lower bounds was within 1 percent of the total energy. This is repeated for as many 
values of Oy and <Py as necessary to evaluate the integrals. 

Unless there is a potential well surrounding the current collection point, there will 
be, for a given value of </? y , a value 6y above which the minimum kinetic energy 
(in eV) equals the potential at that point. Let this value of 0 y be called 9 2 - For low 
plasma temperatures the exponential factor in equation (A9) will change rapidly from 
one for 9y — 9 2 to zero for 6y < 0 2 . The change is less abrupt for higher plasma 
temperatures; nevertheless, Gaussian quadratures tend to give poor results for this 
step-function type of behavior. 

Sufficient accuracy could be obtained by approximating the factor 
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by a function that is zero for 0 V < 0^, rises linearly from Cj at 0j to 
C 2 = (3/2)(i//(r p )/E th ) +1 at 0 V = 0 2 , and is equal to C 2 between 0 2 and tt. This is 
shown in figure 17 where the current collection point is the center of the pinhole and 
where there is no cp dependence. The values of 0^, 0 2 , and Cj depend on <p y and 
E th , and they are obtained by plotting g against 0 V for a given value of <Py. Equa- 
tion (A9) becomes 



This integral is now in a form where Gaussian quadratures will give good results. For 
most of the calculations, however, the current collection point was the center of the pin- 
hole, so the integration over cpy was trivial. 

The next thing to consider is how the particle trajectories are calculated. 

Lagrange’s equations (ref. 12) for a particle moving in a conservative field are 


_d(d] ± \_ 3L _ o 
dt \9qj/ 


(A10) 


where q. is a generalized coordinate (not to be confused with charge on ions) and 

dq. 


% = 


dt 


The Lagrangian for a charged particle in an electromagnetic field is 


L =— mV 2 - q(\p - V • A) 
2 


(All) 


Provision is made for the possibility of a magnetic field, since originally the effect 
of a magnetic field on the current was to be included in the calculations. 

Cylindrical coordinates are used with the disk in the z = 0 plane and the z-axis 
the axis of symmetry. In vector form the velocity V is 
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V = pp + pip <p + zz 


(A12) 


where p, ip, z are unit vectors in the p, cp , and z directions. 
Equation (All) becomes 


L = — (p 2 + p 2 <p 2 + z 2 ) - q[j// - (pA p + piphy + zA z )J 


(A13) 


and (A10) becomes 


3Ap . 3A . 3A Z 

p — - p^3 — - Z 

j ap ap ; 


— (mp + qA ) = mp<p 2 + q<??A - qf-^- - P 

dt p * \3P 

— (mp 2 <p + qpA ) = q (- + P — - + pip — — + z — 

dt * \ dcp d(p d<p d(p/ 

. / dA n dA, n 3A_\ 

— (mz + qA ) = q ( — + p — - + pq? — 21 + z J 

dt z y az az az az f 


r\ 


> (A 14) 




Equations (A14) can be written in terms of the components of the magnetic field B 
by the use of 


B = Vx A 


(A15) 


and 


dA c 3A . 3A . 3A . 3A 
— £ = — 2 + p — 2 + (p — 2 + z — 2 
dt at dp d(p az 


(A16) 


where c stands for any of the components of A. 
In cylindrical coordinates equation (A15) is 


1 aA z 
B = — — - 

p p d<p 

dA 


dA 


<P 


az 


B * = 


P 3A Z 


9z 


dp 


'N 


B 1 3A P 

Z p dp p d(p 


J 


(A17) 
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Using equations (A16) and (A17), equations (A14) become 


mp = q (- + m pep 2 + q 

\ dp at / 


p(p 




3p 




— (mp 2 ^) =q(- 
dt \ 3 ip 


3A 




at 


+ q 


./ aA z 
z p 

W 


3A 


<P 


az 


P 3P . 


P[ A <p + P 


az 


3A 


<P 


3P 


\ az at / 


3A, 


p(^- — 
\ az ap 


\p dep az y 



3<P/ 


mp = m pep 2 + qEp + qificpB ? - zB^) 

j n 

— (mp 99 ) = qpE + qp(zB - pB ) 
dt <P P z 

mz = qE z + q(pB^ - p<?B p ) ' 


(A18(a)) 

(A18(b)) 

(A18(c)) 


Note that E is zero because of symmetry; therefore, in the absence of a magnetic 
field, equation (A18(b)) is d/dt(mp 2 <$?) = 0, which is the statement of conservation of 
angular momentum. If there is no explicit time dependence, E = - V\p, and the equations 
of motion are 


mp = mp<j ? 2 - q-^ + q(p<5?B - zB ) 
3 p z cp 


m cp 






mz = -q^t+ q(pB - pipB ) 

3z V p J 


(A19) 


If the gradient of \p and the components of B are known everywhere, the trajec- 
tory of an incoming ion or election can be reconstructed in the following manner. Let 
the particle land at time t = 0. First, the incoming direction of the particle is re- 
versed so that it will travel backwards over the same path. (It is also necessary to re- 
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verse the direction of B at each point). Next, the position of the particle at a small 
time At earlier is found. If At is small enough, the new coordinates are approxi- 
mately 


P n e* P + P At + -1 p(At) 2 


• 1 2 

(p CX (p+(pb.t+— ip (At) 

n 2 


z z + z At + - z(At) 2 
n 2 _/ 


(A20) 


The initial coordinates p, <p , and z are the coordinates of the current collection point. 
The initial velocities p, ip, and z are calculated from the kinetic energy at impact, 
E ke , and from the impact angles and <p y. The accelerations p, ip, and z are 
found by using equations (A19) and the values of E and B at the current collection 
point. The new velocities are needed before the process can be repeated. These are 


P n ** P + P At') 


^»n ^ ^ + ^ At > 


z„ z + z At J 
n 


(A21) 


The next point is found in the same fashion, with the position and velocity from equa- 
tions (A20) and (A21) used as starting values. The particle is traced step by step back 
along its trajectory to its starting point, either in the space plasma or elsewhere on 
the disk. As a check on the size of At at each step, the kinetic energy at each point is 
calculated from the velocities and compared with the kinetic energy calculated from the 
initial energy of the particle and the potential at that point. The size of At was re- 
duced until the agreement was within 0. 3 percent. This technique has the added advantage 
of making the step size smaller where the field gradients and accelerations are larger. 
This maintains accuracy without wasting computer time. 

The magnetic field can be chosen to represent the field of the Earth or the field of 
some electrical device on the satellite. The only requirement on the magnetic field is 
that it satisfy the relation V • B = 0 to be physically realizable. 

That leaves only the matter of calculating the electric field due to the disk. The 
finite difference approach used by Parker and Whipple was avoided by superposing the 
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potential and electric field due to several disks, each of which was at a uniform voltage 
(fig. 3). 

The voltage around a disk of radius a and uniform voltage is given in oblate 
spheroidal coordinates by Moon and Spencer (ref. 13): 

= — i // 0 cot~*(sinh 77 ) (A22) 


Cylindrical and oblate spheroidal coordinates are related by 


p = a cosh 77 sin 0 ^ 


cp = cp 


z = a sinh 77 cos 9 ) 


(A23) 


The inverse transformations are given by 


cosh 77 = — 
2 a 


y (p + a ) 2 + z 2 + y (p - a ) 2 + 


<p = (p 


(A24) 


sin 6 = — 

2a L 


y(p + a ) 2 + z 2 - y (p - a ) 2 + z 2 


The parameter a enters in the following way: surfaces of constant 77 are oblate 
spheroids; the oblate spheroid for 77 = 0 reduces to a disk of radius a. Since there is 
no built-in inverse contangent function on the computer, it is convenient to rewrite 
equation (A22). Let 7 = (^/2)(i// 1 /i// 0 ) and /3 = sinh 77 so that /3 = cot t which gives 



Substituting for r and /3, 
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1 




1 

cosh 7 ] 


From equation (A24) 

o 

z) =- Sin 


Equation (A25) gives the voltage distribution about a single disk of radius a and voltage 
*o‘ The P rinci Pl e °f superposition is used to construct a solution that represents a 
pinhole surrounded by insulation (fig. 2). To see how this is done, consider the voltage 
at z = 0. From equation (A25) 

*l(P > a, z = 0) = l\f/ 0 sin" 1 ^ (A26) 

V liP s a, z = 0) = (A27) 

For the purpose of illustrating how the solution is obtained, suppose that pinhole- 
insulation system is to have the distribution given by having the three points 

P = r 0 (l)> r 0 (2), r 0 (3) at the voltages V Q (1), V q (2), V q ( 3). From equations (A25) to 
(A27), the general solution is 

3 

X 

Pip, Z)=i ) V. (i) sin" 1 

'L 

i=l 
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where the numbers V b (i) are the superposition voltages and satisfy the equations 

VD = V b (l) + V b (2) + V b (3) I 


V 0 (2) =- V b (l) sin" 1 ^ + V b (2) + V b (3) 
n b r (2) b b 


2 1 r n( 1 ) 2 1 r n( 2 ) 

V 0 (3) = - V b (l) sin" 1 -2_ + £ y ( 2 ) sin’ 1 + V h (3) 
n b r 0 (3) n b r Q (3) b 


(A29) 


The solution given in (A28) is generated by superimposing the voltages of three disks of 
radii r Q (l), r Q (2), and r 0 (3), having, respectively, the constant voltages V b (l), V b (2), 
and V b (3). On the z = 0 plane, i^(p, z=0) automatically reduces to 


*[P <r o (l),0] = V 0 (l)" 
rp[p = r 0 (2), 0] = V 0 (2) l 
V[P = r 0 (3), 0] = V 0 (3) 


(A30) 


The pinhole has radius r Q (l) and voltage V Q (1). The region between r Q (l) and 
r Q (3) is the insulation. In practice, equations (A29) are inverted numerically to get the 
numbers V b , and these are retained in the memory of the computer. The voltage is 
calculated at each point where it is desired by using equation (A28). 

The components of the electric field are still needed. The total electric field is 
calculated by adding the contributions due to individual constant voltage disks, as in 
calculating the potential. This can be seen by using the definition E = -Vjj/, where 4/ 
is given by equation (A28). The expressions for E pl and E z j, the p- and z- 
components of the electric field for a single disk of radius a and voltage V, are most 
easily obtained by converting from oblate spheroidal to cylindrical coordinates. 

The electric field for a single disk expressed in oblate spheroidal coordinates is 



v a 


1 


cosh f] 




cosh 2 r? - sin 2 0 


(A31) 


where rj is the unit vector in the direction of increasing r\. 


Let e„ be defined by 
77 
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(A3 2) 


"V 


e • e 

TJ 7] 


The vector e^ is parallel to rj but is not necessarily a unit vector, 
unit vectors in cylindrical coordinates (ref. 14, p. 120) by 


r =i£p + i5z + i£ * 

■ ” 07} 07] 077 


From equation (A23) 


e^ = a sinh 77 sin 0p + a cosh 77 cos 0z 

— — 2 2 2 2 2 
e^ • e^ = a (sinh 77 sin 0 + cosh 77COS 0) 


= a 2 [(cosh^rj - 1) sin 2 0 + cosh 2 T 7 cos 2 ©] 


a 2 (cosh 2 77 - sin 2 0) 


so that the unit vector 77 is 


77 = 


^/cosh 2 77 - sin 2 0 


(sinh 77 sin 0p + cosh 77 cos 0z) 


From equation (A31), the components of the electric field are 


U _ 2 *0 sinh 77 sin 0 

n a cosh 77(cosh 2 7j - sin 2 ( 


„ 2^o cos 0 

E z " 


7 r a ,2 . 2„ 

cosh 77 - sin 0 


To get an expression in cylindrical coordinates, first use the relation 


It is related to the 


(A3 3) 


(A3 4) 


(A3 6) 
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cosh 2 7 ] - sin 2 e = — \ (p + a) 2 + z 2 \ (p - a) 2 


g _ 2 ^ o a sinh 7) sin 6 
P n cosh tj 


(p + a) 2 + z 2 t/(p - a) 2 + z 2 


2x}/ a. 

E z = cos 6 


(p + a) 2 + z 2 \ (p - a) 2 + z 2 


Equation (A24) supplies the expressions for cosh tj and sin 9, and 


sinh tj = Vcosh tj - 1 


cos 6 = \1 - sin 0 
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APPENDIX B 


^ 1>^2 

E 

E, 


E 


J ke 

ke°° 


E 


mm 


J th 

S zl 

! pi 


V 


j(r p ) 


J eo 


•’10 

k 

L 

m 

m e 

n 

n e 

n o 

P 

q 


SYMBOLS 

magnetic vector potential, V-sec/m 

o 

magnetic field, Wb/m 

parameters used in approximating integral in eq. (A9) 
electric field, V/m 

kinetic energy of particle arriving at current collection point, eV 
kinetic energy of particle in undisturbed plasma, eV 

minimum kinetic energy of particle arriving at current collection point, eV 

thermal energy of particles in undisturbed plasma, eV 

z component of electric field about uniform potential disk 

p component of electric field about uniform potential disk 

vector in 17 direction 

distribution function, (A/m )/(m/sec) 

2 

normal current density at current collection point, A/m 

2 

normal electron current density, A/m 

2 

electron current density in undisturbed plasma. A/m 

2 

normal ion current density, A/m 

2 

ion current density in undisturbed plasma, A/m 

Boltzmann constant, J/K 

Lagrangian function, J 

particle mass, kg 

electron mass, kg 

unit normal at current collection point 

-3 

electron density in undisturbed plasma, m 

-3 

particle density in undisturbed plasma, m 
current collection point 
charge on particle, C 
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- 3 > 


cig charge on electron, C 
th 

Qi i generalized coordinate in Lagrangian function 

Tp position vector of current collection point 

r 0 (i) points at which potential is specified for nonequipotential disk 

T temperature, K 

T g electron temperature, K 

Tj ion temperature, K 

V particle velocity, m/sec 

V|j(i) superposition voltages for nonequipotential disk, V 

V Q (i) specified voltages on nonequipotential disk, V 

Woo particle velocity in undisturbed plasma, m/sec 

r\, 9, ip coordinates in oblate spheroidal coordinate system 

unit normal in tj direction 

^ v an S les that specify direction of incoming particle at current collection point 
(see fig. 15) 

p,<p,z coordinates in cylindrical coordinate system 
p,<p,z unit vectors in p, <p , and z directions 
<MP, z) electrostatic potential, V 
^ p electrostatic potential on pinhole, V 

\p Q uniform potential, V 

^ j(P> z) electrostatic potential due to single disk at Q , V 
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TABLE I. - ION AND ELECTRON CURRENT 


DENSITIES AT VARIOUS LOCATIONS ON 
INSULATOR - VERY POSITIVE CASE 


Thermal 

energy, 

E th> 

eV 

Distance 

along 

radius, 

cm 

Voltage 
at that 
point, 

V 

Electron 

current 

density, 

j e’ 2 

A/m 2 

Ion current 
density, 
h, 

A/m 2 

0.1 

■ 

fl 

375 

1. 18xl0" 2 

0 

1 

■ 

| 



3.74X10' 3 

0 

10 





1.20X10' 3 

2.4X10' 33 

100 

■ 

B 



4. 39x10' 4 

4. 7x10' 8 

0.1 

E 

fl 

225 

7.20X10' 3 

0 

1 


I 



2.28X10' 3 

0 

10 

I 




7.39X10' 4 

8.7X10' 22 

100 

B 

B 



2. 95x10' 4 

6.7xl0' 8 

0.1 

3. 

8 

125 

4. 10X10' 3 

0 

1 





1.30X10' 3 

0 

10 





4.32X10' 4 

3. lxlO' 15 

100 





1.99X10" 4 

3. 0x10' ^ 

0.1 

4. 

8 


2.41X10' 3 

0 

1 



■ 

■ 

7.68X10' 4 

3. lxlO' 57 

10 




1 

2. 62x10' 4 

7. lxlO' 12 

100 



1 

■ 

1.44X10" 4 

6.5X10' 7 

0. 1 

E 

■ 


1.55X10' 3 

0 

1 

B 

■ 


■ 

4.97X10' 4 

9.0X10' 57 

10 

l 

■ 



1.80X10' 4 

1.2x10' 11 

100 

fl 

B 

■ 

1 

1.30X10' 4 

1.9x10" 6 

























TABLE n. - DIMENSIONLESS COLLECTED ELECTRON CURRENT AT 
CENTER OF PINHOLE, WITH INSULATION ADJUSTED TO PRODUCE 
VARIOUS VOLTAGE DISTRIBUTIONS AROUND DISK 
(Disk radius, 5. 1 cm] 


(a) Pinhole radius, 0. 19 cm 


Thermal 


Disk voltage distribution 


energy, 

E th> 

eV 

Very 

Slightly 

Slightly 

Very 

positive 

positive 

negative 

negative 


Collected dimensionless electron current, 

(j e /j eo )(E th/^P^ 

0.1 

1.06 

0.456 

3xl0~ 10 

0 

1 

1.06 

.537 

3xl0' 2 

lxlO" 34 

10 

1.06 

.762 

0.45 

6x10" 5 

100 

1.07 

.974 

.9 

0.3 

1000 

1.27 

1.09 

1. 2 

1.06 


(b) Pinhole radius, 0. 285 cm 


Thermal 


Disk voltage distribution 


energy, 

E th’ 

eV 

Very 

Slightly 

Slightly 

Very 

positive 

positive 

negative 

negative 


Collected dimensionless electron current, 

Meo^W 

0. 1 

1.09 

1.07 

0.749 

0 

1 

1.09 

1.08 

.778 

8xl0" 26 

10 

1.09 

1.08 

.950 

7X10" 4 

100 

1. 11 

1.09 

1.07 

0.4 

1000 

1.23 

1.21 

1. 22 

1. 13 


TABLE HI. - EFFECT OF INSULATION ON PINHOLE CURRENT 


Thermal 

10-kV, 0. 19- cm circular disk 

10-kV, 0. 19- cm circular disk 

energy, 

with no insulation 

plus insulation (very pos. case) 

E th’ 

eV 

CoUected dimensionless electron current, (i e /i eo )( E th /^ p ) 

0.1 

1.41 

1.06 

1 

1.41 

1.06 

10 

1.41 

1.06 

100 

1.42 

1.07 

1000 

1.50 

1.27 




















Thermal energy, E^, eV 

Figure 14 - Collected current density at center pf pinhole. Very positive case, small 
pinhole. 


z 



Figure 15. - Coordinates of current collection point and impact angles of trajectory. 
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